Dynamics of Gauge Fields 



1 



UTCCP-P-84 
UTHEP-421 
April 2000 



Large-Flavor QCD on the Lattice* 



Y. Iwasaki, K. Kanaya, S. Kaya, T. Yoshie 

Institute of Physics (Center for Computational Physics), University of Tsukuba, 
Tsukuba, Ibaraki 305-8577, Japan, iwasaki@rccp.tsukuba.ac.jp, 
kanaya@rccp. tsukuba. ac.jp, kaya@imslab.co.jp, yoshie@rccp. tsukuba. ac.jp 
S. Sakai 

Faculty of Education, Yamagata University, Yamagata 990-8560, Japan. 
sakai@kecszao.kj.yamagata-u.ac.jp 



Abstract 

We study the nature of the QCD vacuum at general number of flavors Np 
by numerical simulations on the lattice. Combining the results with those of the 
perturbation theory, we propose the following picture: 1) For Np > 17, there 
exists only one IR fixed point at vanishing gauge coupling, i.e., the theory in the 
continuum limit is trivial. 2) For 16 > Np > 7, there is a non-trivial fixed point. 
Therefore, the theory is non-trivial with anomalous dimensions, however, without 
quark confinement. 3) For Np < 6, theories satisfy both quark confinement and 
spontaneous chiral symmetry breaking in the continuum limit. 

1. Introduction 



The beta function of QCD is universal up to two-loop in the pertur- 
bation theory: (3(g) = —b g 3 — big 5 , with b = — ^ (ll — an d &i = 
p.~2]2 (l02 — y-^f) ■ When the number of flavors (N F ) exceeds 16|, b becomes 
negative, so that the asymptotic freedom of QCD is lost. In this case, the origin 
g = becomes an IR fixed point, i.e. the low energy limit of QCD is a free 
theory. Quarks cannot be confined and the chiral symmetry cannot be broken 
spontaneously. Np > 16| is a sufficient, but not a necessary condition for these. 
A natural question is, whether confinement and spontaneous breaking of the chiral 
symmetry are satisfied at all Np below 16|. 
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It is well-known that the second coefficient bi changes its sign already at 
Np 8.05. Therefore, a non-trivial IR fixed point may appear at finite g. From 
the two-loop beta function, this happens for Np = 9-16. At least for Np ~ 16, 
the IR fixed point shows up in a perturbative region. Therefore, it is plausible 
that the full beta function has a non-trivial IR fixed point for N' < Np < 16 with 
some N' < 16. When such an IR fixed point exists, the coupling constant cannot 
become arbitrarily large in the IR region — this will imply that quarks are not 
confined. In particular, for Np ~ 16, QCD is perturbative also in the IR limit, 
i.e., we cannot expect confinement. 

In this paper, we study the iVV- dependence of the QCD vacuum. A non- 
perturbative investigation is required. Lattice QCD is the only systematic method 
that enables us to study non-perturbative properties of QCD. We performed a 
series of lattice simulations for a wide range of Np, from 2 to 300, to study the 
phase structure of QCD. When the phase diagram becomes clear, we are able 
to see the nature of the QCD vacuum in the continuum limit, and eventually 
answer the question about the condition on Np for confinement and spontaneous 
breakdown of the chiral symmetry. 

In 1982, Banks and Zaks published a pioneering paper on the iV^-dependence 
of the QCD vacuum [[]]]. Based on an early result obtained on the lattice 0, they 
assumed that, in the strong coupling limit, the theory is confining and the beta 
function is negative, for any value of Np. Using the perturbative results for the 
beta function discussed above, they conjectured Fig. |T](a) as the simplest Np de- 
pendence of the beta function, and studied the phase structure of QCD based on 
this beta function. The assumption of negative beta function in the strong cou- 
pling region leads to an additional non-trivial UV fixed point for Np > N'. Due 
to this UV fixed point, their conjecture for the phase structure is complicated. 
[For different approaches see Refs. || U, H H-] 

In the arguments of Banks and Zaks, the assumption of confinement and 

negative beta function in the strong coupling limit plays an essential role. Here, 

it should be noted that the lattice study on which this assumption was based, is 

a study of a pure gauge theory 0. With dynamical quarks, the vacuum structure 

can be much more complicated. Actually, in the realistic cases, there exist no 

proofs of confinement for general Np even in the strong coupling limit .[] A non- 

1 We can formally rewrite the theory in terms of mesons and baryons. We may argue that, 
when the resulting effective action for hadrons is a well-defined action of weakly interacting 
particles, it is plausible that quarks are confined — the hadrons are good variables to describe 
the dynamical contents of the theory. In the strong coupling limit, the effective action can be 
computed either in the large N c limit 0, using a meanfield approximation (1/d expansion) 
H D, or using a heavy quark mass expansion IJ, |T^]. In these cases, the effective action 
for mesons and baryons seems to lead to a well-defined world with spontaneously broken 
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b) 

Fig. 1. Renormalization group beta function of QCD. (a) Conjecture by Banks 
and Zaks [0J assuming confinement in the strong coupling limit for all Np. (b) 
Our conjecture deduced from the results of lattice simulations. 

perturbative investigation on the lattice is required. 

In a previous study, we performed simulations of QCD in the strong cou- 



pling limit at various Np 111 1|| ■ We found that, when Np > 7, light quarks are 
not confined and chiral symmetry is restored, even in the strong coupling limit. 
Here, we extend this study to weaker couplings and to a wider range of N F up 
to 300. Combining the lattice results and perturbative arguments, we conjec- 
ture Fig. 0(b) for the non-perturbative beta function of QCD. This leads to the 
following ^-dependence of the QCD vacuum. 

• When Np is smaller than a critical value, the beta function is negative for all 
values of g. Quarks are confined and the chiral symmetry is spontaneously 
broken at zero temperature. We find that the critical number of flavors is 
6. (Corresponding critical value of Np is 2 for N c = 2, to be compared with 
the two- loop value 5 [I2"|.) 



chiral symmetry, suggesting quark confinement. For the realistic cases of N c = 3 with light 
quarks, however, the effective action is quite complicated; diagrams with quark loops, which 
are suppressed at large N c , large d or large m q , become important. Therefore, in the realistic 
cases, it is difficult to analytically deduce the vacuum structure as a function of Np even in 
the strong coupling limit. 



4 



On the other hand, when Np is equal or larger than 17, we conjecture that 
the beta function is positive for all g, in contrast to the conjecture by Banks 
and Zaks shown in Fig. [L](a). The theory is trivial in this case. 

When Np is between 7 and 16, the beta function changes sign from negative 
to positive with increasing g. Therefore, the theory has a non-trivial IR 
fixed point, i.e., the theory is non-trivial with anomalous dimensions. We 
conjecture that quarks are not confined in this case. 



Our results are in part reported in Refs. ]T2[ and JIB] . 

This report is organized as follows: In Sec. |2], we describe our lattice model 
and simulation parameters. Results for the phase structure in the strong coupling 
limit are given in Sec. |3j[ The phase structure at finite coupling is discussed in 
Sec. [|]. Finally, in Sec. [5J, we summarize our conclusions and conjectures. 

2. Model and simulation parameters 

We consider a 4-dimensional Euclidian lattice with the lattice spacing a 
and finite lattice volume (N s a) 3 x (N t a), where N s and N t are the number of 
lattice sites in spatial and temporal directions. The action consists of the gauge 
part and the quark part, S = S gauge + S quar k- For the gauge part, we use the 
standard one plaquette action. 



Sgauge — „ E Tr (Ux^Ux+fi^U^^fi^U^^Jj , (1) 



where U xli is the gauge connection between the site x and the neighboring site in 
the positive \i direction, and (3 = 6/g 2 . 

For the quark part, there are several alternatives. Two conventional choices 
are the staggered (Kogut-Susskind) fermion action |14j and the Wilson fermion 



action fL0| . Staggered quark has been preferably used in previous studies |15 
because a part of the chiral symmetry is preserved in this formulation. On the 
other hand, the flavor structure is quite complicated: The action is local only 
when Np is a multiple of 4. Off the continuum limit, there exists a lattice artifact 
mixing these 4 degenerate flavors. For the more realistic cases Np = 2 and 3, we 
modify by hand the power of the fermionic determinant in the numerical path- 
integration. This necessarily makes the action non-local, that sometimes poses 
conceptually and technically difficult problems. 

Another conventional choice for the quark action is the Wilson fermion 
action, which we adopt in the followings. 

S quark = E - K {$ f x {l-l»)U x ,^l +[l + ^{l+lM,M}] > ( 2 ) 

f,x 
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The index / (= 1, • • • , Np) is for the flavors. We assume that all Np quarks are 
degenerate with the bare mass mo, which is related to the hopping parameter K 
by niQ = (l/2a)(l/ K — 1/K C ), where K c = 1/8 is the point where the bare mass 
vanishes. In this formulation, the chiral symmetry is explicitly broken at finite a. 
On the other hand, the flavor symmetry as well as the C, P and T symmetries 
are manifestly satisfied without mixing etc. also on the lattice. 

In this study, we use the Wilson fermion formalism because this is the 
only formalism that preserves manifest flavor symmetry. This feature will be 
important in a study of Np- dependence in QCD. Concerning the chiral symmetry, 
a perturbative study of Ward identities shows that the effects of the 0(a) chiral 
breaking terms can be absorbed by appropriate renormalizations, including an 



additive renormalization of quark mass fL6 |. As a result, the location of the chiral 



limit K c shifts from the tree value 1/8 as a function of (3. 

In order to take the additive mass renormalization into account, we define 
the current quark mass m q in terms of an axial vector Ward identity Jl7, [16 . 



2m,(0|P|7r(p = 0)) = -ffl, ( | A 4 \ n(p = 0) ) (3) 

where P is the pseudoscalar density, A4 the fourth component of the local axial 
vector current, and m n the pion (screening) mass. A multiplicative renormaliza- 
tion factor for the axial current, which is not important in this study, is absorbed 
into the definition of the quark mass. We then define K c by the condition m q = 0. 
Numerical simulations for small Np shows that K c (f3) is a smooth curve connect- 
ing 1/8 at p = 00 and 1/4 at = 0. 

In QCD, even when the confinement is realized at low temperatures, we 



expect deconfinement at high temperatures []18| . Therefore, we have to study 
the temperature dependence of the results to distinguish between the finite tem- 
perature deconfinement transition/crossover due to high temperatures and a bulk 
transition at zero temperature due to the effects of many flavors. On a lattice with 
N t sites in the Euclidian time direction, the temperature is given by T = 1/N t a. 
When the beta function is negative, as in the case of small Np, the lattice spacing 
a is a decreasing function of (3 = Q/g 2 . In this case, T increases with increasing 
(3 for a fixed finite Np On the other hand, when the beta function is positive, 
T decreases with increasing f3. Therefore, in order to study the temperature 
dependence for general Np, we have to compare the results at various N t . 

Numerical study shows that the value of m q defined through the axial 
vector Ward identity does not depend on whether the system is in the high or 
the low temperature phase when (3 is large; /3 > 5.5 for Np = 2 |T9| . When we 
define the pion decay constant f n by ( | A4, \ Tr(p = 0) ) = m^f^, either m n = or 
f n = has to be satisfied in the chiral limit. The chiral symmetry is restored when 
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Fig. 2. Plaquette at /3 = as a function of 1/K. (a) For jVp = 7-300 at N t = 4. 
(b) For N F = 240 at N t = 4-16. 



f w = in the chiral limit, while = if the chiral symmetry is spontaneously 
broken [2(J . 



We perform simulations on lattices 8 2 x 10 x N t (N t = 4, 6, and 8), 
16 2 x 24 x N t (N t = 16) and 18 2 x 24 x N t (N t = 18). We vary N F from 2 to 
300. For each Np, we study the phase structure in the coupling parameter space 
((5, K). It should be noted that, in QCD with dynamical quarks, there are no order 
parameters for quark confinement. We discuss about confinement by comparing 
the screening pion mass and the lowest Matsubara frequency, and, simultaneously, 
consulting the values of plaquette and the Polyakov loop. In the followings, we 
call the pion screening mass simply the pion mass, and similarly for the quark 
mass. Further details about our simulation are described in Refs. fi~2 , [13| 



3. Strong coupling limit (5 = 



In a previous paper we have shown that, when Np > 7, light quarks 
are deconfined and chiral symmetry is restored at zero temperature even in the 
strong coupling limit. For iV^ < 6, we have only one confining phase from the 
heavy quark limit K = up to the chiral limit K c = 0.25. On the other hand, 
for Np > 7, we find a strong first order transition at K — Kd- 

• When quarks are heavy (K < Kd), both plaquette and the Polyakov loop 
are small, and m n satisfies the PCAC relation m 2 oc m q . Therefore, quarks 
are confined and the chiral symmetry is spontaneously broken in this phase. 
We find that m q in the confining phase is non-zero at the transition point 
Kd- While an extrapolation in 1/K shows that m 2 decreases towards zero 
at K = 0.25, the chiral limit does not belong to this phase (Kd < K c ). 
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Fig. 3. The transition point at (3 — versus Np for N t = 4 and N t > 8. 

For clarity, data at Nt — 8 for N F = 300 is slightly shifted to a larger Np in 
the figure. 



On the other hand, when quarks are light (K > Ka), plaquette and the 
Polyakov loop are large. In this phase, m n remains large in the chiral limit 
and is almost equal to twice the lowest Matsubara frequency ir/N t . This 
implies that the pion state is a free two-quark state and, therefore, quarks are 
not confined in this phase. The pion mass is nearly equal to the scalar meson 
mass, and the rho meson mass to the axial vector meson mass. The chiral 
symmetry is also manifest within corrections due to finite lattice spacing. 



See 11 for details. 



We extend the study to larger Np. We show the results of plaquette at 
Nt — 4 for Np = 7-300 in Fig. 0(a). Clear first order transition can be seen at K 
below 0.25 (1/K > 4). We then study the N t dependence of the results, as shown 
in Fig. 0(b) for the case of Np = 240. We find that, although the transition point 
shows a slight shift to smaller 1/K when we decrease T = 1 / N t a from N t = 4 to 
8, it stays at the same point for N t > 8. The chiral limit, 1/K = 4, cannot be 
achieved in the confining phase even at T = 0. Similar result for Np = 7 is given 
in Ref. fTI]] . Therefore, we conclude that the transition is a bulk transition. 

From these results, we obtain the phase diagram at (3 — shown in Fig. 0. 
The first-order transition line separating the confining phase and the deconfining 
phase reaches 1/K — 4 (chiral limit for small Np) between Np = 6 and 7. Note 
that the critical value of the inverse critical hopping parameter is greatly increased 
from 1/K = 4.08 to 7.75 when Np increases from 7 to 240. Recall that 1/K = 8 
is the point for massless free quarks at (3 = oo. 
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Fig. 4. The phase structure for Np < 6 for (a) T > 0, and (b) T = 0. 




Fig. 5. Phase diagram for (a) Np = 12 and (b) Np = 18. Dark shaded lines 
represent our conjecture for the bulk transition line in the limit N t = oo. 



4. Phase structure at finite (3 

Let us now study the case f3 > 0. We perform simulations at (3 — 0.0-6.0 
varying Np from 2 to 300. On a lattice with a fixed finite Nt, we have the finite 
temperature deconfining transition when the quark mass is sufficiently heavy. In 
particular, at K = 0, we have the first order deconfinement transition of the pure 
gauge theory at finite (3 irrespective to the value of ./V^. 

For Np < 6, the finite temperature transition/crossover line crosses the 
K c line at finite (3 |[20|| . A schematic diagram of the phase structure for this case 
is shown in Fig. [I]. (In the followings, we call the transition/crossover line simply 
as transition line.) We have a monotonous RG flow towards the strong coupling 
limit also on the K c line. Therefore, the whole finite temperature transition line 
moves towards a larger (3 as N t is increased. In the limit N t = oo (i.e., at T = 0), 
only the confining phase is left. 

For Np > 7, we have seen that there exist a bulk deconfining transition 
in the strong coupling limit. We have to clarify the relation between the finite 
temperature transition at small K and the bulk transition at (3 — 0. From our 
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Fig. 6. Phase diagram for Np = 240. 



simulations, we obtain the phase diagrams sown in Fig. [5] for Np = 12 and 
18, and in Fig. for Np = 240. When N t = 4 or 8, the phase boundary line 
between confining and deconfining phase bends down at finite (3 due to the finite 
temperature phase transition of the confining phase. The darker shaded lines 
are our conjecture for the phase boundary between the confining phase and the 
deconfining phase at zero temperature (Nt = oo). 

In Figs. |5] and the dashed lines in the deconfining phase represents the 
chiral limit, m q = 0. This line also corresponds to the minimum point of m^. We 
have checked that, for the case Np = 240, the quark propagator in the Landau 
gauge actually shows the chiral symmetry, 75^(2)75 = —G(z), on the m q = 
line ||12|| . For Np> 240, we find that the m q = line in the deconfining phase, 



which starts at 1/K — 8 in the weak coupling limit (3 = 00, reaches the strong 
coupling limit, as shown in Fig. |6] for Np = 240. For a smaller Np < 100, because 
the bulk transition line in the strong coupling region shifts toward larger K 
with decreasing Np, the m q = line in the deconfining phase hits the line at 
finite f3. For example, in the case Np = 18, it hits at /3 — 4.0 - 4.5, as shown in 
Fig. 0(b). 



4.1. RG-flow for N F > 17 

The m q = point at (3 = 00 is a trivial IR fixed point for iV^ > 17. The 
phase diagram shown in Fig. [T] suggests that there are no other fixed points on 
the m q = line at finite f3. 

Figure 0(a) shows the results of m% and 2m q in the deconfining phase for 
Np = 240 at Nt = 4. We find that the shape of and 2m q as a function of 
1/K only slightly changes for 1/K < 8 when the value of (3 decreases from 00 
down 0. We obtain similar results also for N t = 8. This suggests that quarks are 
almost free down to (3 = in the deconfining phase. The results for iV^ = 300 
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are essentially the same, except for very small shifts of the transition point. 

We make a Monte Carlo renormalization group (MCRG) study along the 
m q = at Np = 240. Performing a block transformation with scale factor 2, we 
estimate the quantity A/3 = (3 (2a) — (3(a): We generate configurations on an 8 4 
lattice on the m q = points at (3 — and 6.0 and make twice blockings. We 
also generate configurations on a 4 4 lattice and make once a blocking. Then we 
calculate A/3 by matching the value of the plaquette at each step. We obtain 
Af3 ~ 6.5 at (3 — and 10.5 at (3 — 6.0.0 The value obtained from the two-loop 
perturbation theory is Af3 ~ 8.8 at (3 = 6.0. The signs are the same and the 
magnitudes are comparable. This suggests that the direction of the RG flow on 
the m q = line at (3 = and 6.0 is the same as that at (3 = oo. This further 
suggests that there are no fixed points at finite (3. These imply that the theory is 
trivial for Np = 240. 

The area of the deconfining phase decreases with decreasing Np in the 

strong coupling region. However, the shape and the values of m% and m q are 

2 In order to get a more precise value of A/3, one has to make many steps of block trans- 
formations, together with a fine tuning of parameters of the block transformation, and do 
matching using several types of Wilson loops. We reserve elaboration of this point for future 
works. For Np = 240, because the velocity of the RG flow is large, we will be able to obtain 
the correct sign and approximate value of A/3 by a simple matching. 
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quite similar when we vary Np from 300 down to 17. Fig. [T](b) shows and 
2m q for N F = 18 - 300 at /3 = 6.0 and 0. At (3 = 6.0, the shapes of m\ are 
almost identical to each other, except for a small shift toward smaller 1/K as Np 
is decreased. 

These facts suggest that, for Np > 17, the nature of physical quantities 
in the deconfining phase is almost identical to that observed for Np = 240 and 
300, i.e., quarks are almost free. Therefore, we conjecture that the direction 
of the RG flow at (3 <Q is identical to the case Np = 240. Combining this 
with the perturbative result that (3 = oo is an IR fixed point for Np > 17, we 
conjecture that the RG flow along the massless quark line in the deconfining phase 
is uniformly directing to smaller (3, as in the case of Np = 240. When this is the 
case, the theory has only one IR fixed point at f3 — oo, i.e. the theory is trivial 
for Np > 17. 

4.2. 16 > Np > 7 

As discussed in Sec. [T], quark confinement is lost for Np > 7 at (3 = 0. 
We have performed simulations for Np = 12 and 7. The phase diagram for 
Np = 12 is shown in Fig. |5](a). At /? < 6.0 we studied, the gross feature of the 
phase diagrams for Np = 7 and 12 is quite similar to the case Np = 18 shown in 
Fig. 0(b). Physical quantities at (3 < 6.0 also show behavior similar to that shown 
in Fig. [7] for Np > 17. Therefore, we consider it probable that the direction of the 
RG flow in the deconfining phase at small (3 is towards a larger (3 as in the case 
of Np > 17. However, the direction of the RG flow at (3 = oo is opposite to that 
for Np > 17 because the theory is asymptotically free for Np < 16. This means 
that we have an IR fixed point somewhere at a finite value of (3. The continuum 
limit is governed by this IR fixed point. 

On the m q = line around (3 = 4.5, we find that m n is roughly twice the 
lowest Matsubara frequency. This implies that the quarks are not confined and 
almost free. Therefore, the anomalous dimensions are small at the IR fixed point, 
suggesting that the IR fixed point locates at finite (3. Unfortunately, as shown 
in Fig. [T], the massless line hits the boundary in the coupling parameter space. 
Without much space in the strong coupling region where simulation is easy, it is 
hard to study numerically the location of the IR fixed point. A detailed study of 
RG flow in a wider coupling parameter space is required. We reserve these studies 
for future works. 
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Fig. 8. Phase structure and the RG flow at T = for general number of flavors. 
Thin dashed lines in the phase diagrams for Np > 7 represent the location of 
the chiral limit for Np < 6 guide for eyes. 



5. Conclusions 

Summarizing our numerical results and conjectures, we propose Fig. [8] for 
the phase structure of QCD at T = 0: When N F < 6, quarks are confined for 
any values of the current quark mass and f3. The chiral limit K c (f3), where the 
current quark mass m q vanishes, belongs to the confining phase. When Np > 7, 
there is no chiral limit in the confining phase. There is a line of a first order phase 
transition from the confining phase to a deconfining phase at a finite current quark 
mass for all values of (3. The chiral limit exists only in the deconfining phase, and, 
therefore, can not be taken from the confining phase. 

In order to see the orientation of the RG flow on the chiral limit line in 
the deconfining phase, we performed a MCRG study at Np = 240, because, in 
this case, we have sufficiently long massless quark line in the strong coupling 
region where numerical simulation is easy. We found that the RG flow along the 
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massless quark line at small (3 is consistent with the result of the perturbation 
theory. This suggests that the direction of the RG flow is uniform for all values 
of /3, as shown in Fig. R], Therefore, there is only one IR fixed point at (3 = oo, 
so that the theory in the continuum limit is trivial. Accordingly, we find that, 
in the deconfining phase, the i^-dependences of pion screening mass and current 
quark mass at small (3 are almost identical to those of free Wilson quarks. 

Also for smaller values of Np > 7, the general features of the phase diagram 
and the independence of physical quantities are quite similar to the case of large 
Np ~ 240. Therefore, we conjecture that the direction of the RG flow is the same 
as in the case of Np = 240, at least for (3 < 6.0 we have studied. 

Together with the results of the perturbation theory, we conjecture the 
following picture. 

• For Np > 17, because the RG flow in the weak coupling limit has the same 
direction due to the lack of asymptotic freedom, we conjecture a uniform RG 
flow towards the trivial IR fixed point, as in the case of Np = 240. Thus the 
theory is trivial for Np > 17. 

• When 7 < Np < 16, because the theory is asymptotically free, the direction 
of the RG flow is opposite at large (3. Therefore, we expect a non-trivial 
IR fixed point at finite (3. In this case, the theory in the continuum limit 
is a non-trivial theory with anomalous dimensions. Because the coupling 
constant in the IR limit is finite, quarks are not confined in this case. 

• When Np < 6, quarks are confined and the chiral symmetry is spontaneously 
broken for any values of the current quark mass and (3. The chiral limit 
belongs to the confining phase. 
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